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Chapter 1 

Feedback Processes at Cosmic Dawn 

Ferrara, A. & Salvaterra, R. 

SISSA/ISAS, Via Beirut 4, 34OI4 Trieste, Italy 



The word "feedback" is by far one of the most used ones in modern cos- 
mology where it is appHed to a vast range of situations and astrophysical 
objects. However, for the same reason, its meaning in the context is often 
unclear or fuzzy. Hence a review on feedback should start from setting the 
definition of feedback on a solid basis. We have found quite useful to this 
aim to go back to the Oxford Dictionary from where we take the following 
definition: 

Fee'dback n. 1. (Electr.) Return of fraction of output signal from 
one stage of circuit, amplifier, etc. to input of same or preceding stage 
(positive, negative, tending to increase, decrease the amplification, etc). 
2. (Biol., Psych., etc) Modification or control of a process or system by its 
results or effects, esp. by difference between the desired and actual results. 

In spite of the broad description, we find this definition quite appro- 
priate in many ways. First, the definition outlines the fact that the concept 
of feedback invokes a back reaction of a process on itself or on the causes 
that have produced it. Secondly, the character of feedback can be either 
negative or positive. Finally, and most importantly, the idea of feedback is 
intimately linked to the possibility that a system can become self-regulated. 
Although some types of feedback processes are disruptive, the most impor- 
tant ones in astrophysics are probably those that are able to drive the 
systems towards a steady state of some sort. To exemplify, think of a 
galaxy which is witnessing a burst of star formation. The occurrence of 
the first supernovae will evacuate/heat the gas thus suppressing the star 
formation activity. Such feedback is then acting back on the energy source 
(star formation); it is of a negative type, and it could regulate the star 
formation activity in such a way that only a sustainable amount of stars is 
formed (regulation) . However, feedback can fail to produce such regulation 

1 
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either in small galaxies where the gas can be ejected by the first SNe or 
in cases when the star formation timescale is too short compared to the 
feedback one. As we will see there are at least three types of feedback, and 
even the stellar feedback described in the example above is part of a larger 
class of feedback phenomena related to the energy deposition of massive 
stars. We then start by briefly describing the key physical ingredients of 
feedback processes in cosmology. 

1.1 Shock Waves 

A shock represents a "hydrodynamic surprise" , in the sense that such a 
wave travels faster than signals in the fluid, which are bound to propagate 
at the gas sound speed, Cg. It is common to define the shock speed, ^,5, in 
terms of c^ by introducing the shock Mach number, M: 

v,^Mcs = M{-fP/pY/\ (1.1) 

where P and p are the gas pressure and density, respectively. The change 
of the fluid velocity v is governed by the momentum equation 

dv 11 

p—+pvVv = -VP-—VB'^ + —B-VB, (1.2) 

ot Stt 47r 

where B is the magnetic field. The density is determined by the continuity 
equation 

%+vVp^-p\/v. (1.3) 

ot 

1.1.1 Hydrodynamics of Shock Waves 

The large-scale properties of a shock in a perfect gas, with B ~Q are fully 
described by the three ratios W2/W1 P2I p\ and P2/ Pi^ determined in terms 
of the conditions ahead of the shock (i.e. vi, pi and Pi) by the Rankine- 
Hugoniot jump conditions resulting from matter, momentum and energy 
conservation: 

pivi ^ P2V2, (1.4) 

Pi+Pivl = P2 + P2vl, (1.5) 

«2(x/52W2 + U2) - Vl{-pivl + Ui) = ViPi - V2P2, (1.6) 

where U is the internal energy density of the fluid. If the fluid behaves as 
a perfect gas on each side of the front we have U = P/{-f — 1). 

In the adiabatic limit it is sufficient to consider the mass and momen- 
tum equations only. By solving such system we obtain: 

V2 Pi 7^1, 2 1 



vi p2 7 + 1 7 + 1M2' 



(1.7) 
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where the Mach number is defined as JVP — v'l/c\ and ci is the sound 
velocity in the unperturbed medium ahead of the shock. For strong shocks, 
i.e. M > 1 and 7 = 5/3, P2IP1 = 4 and 

T. = |i»^ (1.8, 

where Vs — vi — V2 

In the isothermal limit, i.e. when the cooling time of the postshock 
gas becomes extremely short and the initial temperature is promptly re- 
established, 

ei = d^M^ (1.9) 

where Cg is the sound velocity on both sides of the shock. While in the 
adiabatic case the compression factor is limited to four, in an isothermal 
shock much larger compressions are possible due to the softer equation of 
state characterizing the fluid. 

1.1.2 Hydromagnetic Shock Waves 

Assuming B is parallel to the shock front (so that B ■ VB = 0), eq. (|1.5|l 
can be replaced by 

P,, „„; + (!). ft + ,,„5 + (|i), (1,10) 

Magnetic flux conservation requires i?i/pi = B2/P2, so that in the 
adiabatic limit the magnetic pressure B'^ /Sir increases by only a factor 16. 

In the isothermal limit the compression depends only linearly on the 
Alfvenic Mach number Ma 

Pl^^J!^^Ma, (1.11) 

Pi Va,l 



where Va = B/y/Anpi is the Alfven velocity. Hence, some of the shock 
kinetic energy is stored in the magnetic field lines rather than being used 
to compress the gas as in the field-free case; in other words a magnetized 
gas is less compressible than an unmagnetized one. 

1.1.2.1 Structure of Radiative Shocks 

The structure of a shock can be divided into four regions (McKee & Draine 
1991). The radiative precursor is the region upstream of the shock in which 
radiation emitted by the shock acts as a precursor of the shock arrival. The 
shock front is the region in which the relative kinetic energy difference of 
the shocked and unshocked gas is dissipated. If the dissipation is due to 
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collisions among the atoms or molecules of the gas, the shock is collisional. 
On the other hand, if the density is sufficiently low that collisions are 
unimportant and the dissipation is due to the collective interactions of the 
particles with the turbulent electromagnetic fields, the shock is collisionless. 
Next comes the radiative zone, in which collisional processes cause the gas 
to radiate. The gas cools and the density increases. Finally, if the shock 
lasts long enough, a thermalization region is produced, in which radiation 
from the radiative zone is absorbed and re-radiated. 

1.1.3 Supernova Explosions 

In the explosion of a supernova three different stages in the expansion may 
be distinguished. In the initial phase the interstellar material has little 
effect because of its low moment of inertia; the velocity of expansion of 
the supernova envelope will then remain nearly constant with time. This 
phase terminates when the mass of the swept up gas is about equal to the 
initial mass Me expelled by the supernova, i.e. {4TT/3)r^pi = Me, where 
pi is the density of the gas in front of the shock and r^ is the radius of 
the shock front. For Me = 0.25 M© and pi = 2 x lO"^'^ g/cm^, r^ ~ 1 
pc, which will occurs about 60 yr after the explosion. During the second 
phase (Sedov- Taylor phase) the mass behind the shock is determined pri- 
marily by the amount of interstellar gas swept up, but the energy of this 
gas will remain constant. When radiative cooling becomes important (ra- 
diative phase), the temperature of the gas will fall to a relatively low value. 
The motion of the shock is supported by the momentum of the outward 
moving gas and the shock may be regarded as isothermal. The velocity of 
the shell can be computed from the condition of momentum conservation 
(snowplow model). This phase continues until Vg = Taax{vt,Cs}, where vt 
is the turbulent velocity of the gas, when the shell loses its identity due to 
gas random motions. 

1.1.3.1 Blastwave Evolution 

To describe the evolution of a blastwave, the so-called thin shell approxi- 
mation is frequently used. In this approximation, one assumes that most 
of the mass of the material swept up during the expansion is collected in a 
thin shell. We start by applying (we neglect the ambient medium density 
and shell self-gravity) the virial theorem to the shell: 



2 dt 



2 - 2E^2{Ek + Et), (1.12) 



where / is the moment of inertia, and E^ {Et) is the kinetic (thermal) 
energy of the shell. Let us introduce the structure parameter K in the 
following form: 



1 '■« 
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K = —TFr- / dmrv, (1.13) 

MRV2J0 ^ ' 

where M is the ejected mass, R the shell radius. Substituting in equa- 
tion (|1.12|) with V = V2{r/R), we obtain 

Kj^iMrv,) = 2Eiv,/v2). (1.14) 

To recover the Sedov- Taylor adiabatic phase we integrate up to R and 
assume p{R) = po-R~", E = const., K ^ {6 - 2m)/(7 - 2m) 



i?(t) = 



{5-m)E 



-1 l/(5-m) 



^2/(5""). (1.15) 



npoK 

The previous expression can be obtained also via a dimensional approach. 
In fact, in the adiabatic phase the energy is conserved and hence 

E ex A/w^ (X poR^-"'t-^, (1.16) 

so that 

Po 
the fraction of the total energy transformed in kinetic form is 



1,,,, /3 2\ 3E 



(LIE 



2 V4 V 2(5 -m)' 

We note that in the limit of a homogeneous ambient medium {m = 0) we 
recover the standard Sedov- Taylor evolution R (x t^/^ , with 30% of the 
explosion energy in kinetic form. 

A similar reasoning can be followed to recover the radiative phase 
behavior. In this phase V2 = Vs, i.e. the shock wave is almost isothermal 
and K — 1, i.e. the shell is thin. Under these assumptions the energy is 
lost at a rate 

E = -47ri?2 Qpo«?) • (1.19) 

There are two possible solutions of equation (|1.12|) in this case: (a) 
the pressure-driven snowplow in which 

i?(i)cxt2/(7-™)^ (1.20) 

or the momentum-conserving solution 

i?(t)cxti/(4-'"). (1.21) 
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In both cases, the expansion proceeds at a lower rate with respect to the 
adiabatic phase. 

In order to study the evolution of an explosion in a more realistic 
way it is necessary to resort to numerical simulations as the one shown 
in Figure 11.11 From there the inner structure of the remnant is clearly 
visible: indeed most of the mass is concentrated in a thin shell (located 
at about 80 pc after 0.5 Myr from the explosion); a reverse shock is also 
seen which is proceeding towards the center and thermalizing the ejecta up 
to temperatures of the order of 10^ K. In between the two shocks is the 
contact discontinuity between the ejecta and the ambient medium, through 
which the pressure remains approximately constant. 
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Figure 1.1. The structure of spherical supernova remnant calculated us- 
ing a high-resolution numerical hydrodynamical code. The remnant age is 
t = 5 X 10^ yr. This Figure is taken from Cioffi, McKee & Bertschiger (1988) 



1.2 Photoionization 

The ionization balance equation around a source is given by 
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.2dR , 4_ ,_,^ 



AnR^n— + -nR^n^a = 47rS',(0), (1.22) 

where S'i(O) is the ionizing photon rate from the source, n is the neutral 
gas density and a is the recombination coefficient. The solution of the 
ionization balance equation is 

R{t) = Rs[l - exp(-2nat/3)]i/^ (1.23) 

As time goes to infinity the radius tends to the Stromgren radius 

\ an 
The above equations implies that Rg is reached in 



^^-(^^E^^ ■ (1-24) 



ts = -na ~ 2 X lO^n yr, (1.25) 

and the typical ionization front (I- front) velocity is vj ex Si{0)l R?n. 

Depending on the value of vi we can distinguish three types of I-fronts, 
where Ci (c„) is the sound speed in the ionized (neutral) gas: 

(i) R(arefied)-TYPE: u/ > (q,c„) 

(u) Strong D(ense)-TYPE: c„ < w/ < a 

(iii) Weak D(ense)-TYPE: vi < (q,c„) 

In the case of R-TYPE I-front there is no hydrodynamical reaction to 
the passage of the ionizing front, since the its motion is supersonic relative 
to the gas behind as well as ahead of the front. In the case of D-TYPE 
fronts the dynamical timescale becomes short enough for the gas to react 
to photoionization and a shock will form and precede the ionizing front 
thus increasing the ambient gas density. 

1.2.0.2 Photoionization Equilibrium 

It is convenient to divide the ionizing flux into two parts, a 'stellar' part 
S{iy), resulting directly from the input radiation from the star, and a 'dif- 
fuse' part Sdiv), resulting from the emission of the ionized gas, that in a 
pure H gas is due to recombinations into 1^ S state. 

St{y)^S{p) + Sd{p). (1.26) 

Both the recombination radiation and photoionization cross section (ct cx 
v^) are peaked at the threshold, v^. We can identify two approximate cases 
(Ostcrbrok 1989): 
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(i) MENZEL CASE A: optically thin nebula, i.e. t{i^ = i^l) = nH(y{i'L)R < 
1, and the diffuse photons can escape; 

(ii) MENZEL CASE B: optically thick nebula, i.e. diffuse photons are 
absorbed very close to the point at which they are generated ('on-the- 
spot' approximation). 

Menzel case B is usually a good approximation for many astrophysi- 
cal purposes because the diffuse radiation-field photons have t^ ~ i^^, and 
therefore a short mean free path before absorption. Under the hypothesis 
of the on-the-spot approximation, the recombination coefficient is given by 

aeff = a- ais, (1-27) 

resulting roughly in a 60% decrease. The physical meaning is that in opti- 
cally thick nebulae, the ionizations caused by stellar radiation-field photons 
are balanced by recombinations to excited levels of H, while recombinations 
to the ground level generate ionizing photons that are absorbed elsewhere 
in the nebula but have no net effect on the overall ionization balance. 

1.3 Thermal Instability 

A final crucial ingredient to understand feedback and self-regulation is rep- 
resented by the thermal instability. Let us rewrite the usual mass, momen- 
tum and energy equations in a slightly different form (Balbus 1995): 



dlnp 
dt 

' dt 



pV-v^O, (1.28) 

P^ = -VP, (1.29) 



3 dlnP r/o , , , , 

^P^^p-"^^-pL{p,T), (1.30) 

where L is the net cooling function (i.e. heating-cooling) and is taken to be 
a function of density and temperature; in the equilibrium state v = L = Q. 
We perturb the equilibrium by introducing small disturbance. We write 
such perturbations as 5q ex Qif^^+nt).^ ^^ jg then the growth rate of the 
instability. The linearized system yields the dispersion relation 

„3 + ^reT,pn^ + k^cln + \k'^clTQT,p = 0, (1.31) 

where 6 = {plP)L and Qa,b = {dQ/dA)B. 

For large k (small wavelengths) the instability occurs isobarically, that 
is pressure is kept approximately constant as the density increases following 
the process of condensation and cooling. In this case such equations has 
three solutions: 
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nx = -\tQt.p, (1.32) 

5 

2 

712,3 = ±ikcs - —TQt.s- (1.33) 

15 

It is straightforward to recognize that ri2,3 are the modified sound wave 
modes, whereas ni corresponds to the condensation mode. Note also that 
the sign of the temperature derivative of the cooling function determines 
whether or not a mode is unstable (instability occurs if Qt.p < 0). 

In the opposite (isochoric) limit characterized by small k (large wave- 
lengths), the condensation mode has the following expression: 

ni--^reT,p, (1.34) 

Although formally a solution of the dispersion equation, the thermal in- 
stability cannot physically occur under isochoric conditions. The reason is 
that on such large scales there is not enough time to establish pressure equi- 
librium as the cooling time is much shorter than the sound crossing time 
over a wavelength. Thus the density remains approximately constant while 
the gas cools, and pressure equilibrium is likely to be restored promptly 
and rather dramatically by a propagating shock front. 

To better understand the nature of the thermal instability let us plot 
the curve of the equilibrium solution L = in the P — T plane for a hypo- 
thetical cooling function (Figure [TT^ . Above the curve, heating dominates 
over cooling; below L = 0, the gas is cooled. Between Pmin and Pmax a 
constant pressure line will intersect the curve in three points (A, B, and 
C); A and B are thermally stable. This is readily understood since a small 
isobaric increase (decrease) in temperature at either of these points takes 
the system into the cooling (heating) region of the P ~T plane from where 
it is immediately driven back to equilibrium position. On the other hand, a 
small increase (decrease) in the temperature in the unstable point B takes 
one into the heating (cooling) region, with the result that the system is 
driven yet further from its equilibrium. This leads to the coexistence of 
two (or more) thermodynamic phases and provides the main argument to 
interpret the multiphase medium commonly observed in galaxies and in 
cosmic gas. We will see later on how this concept is closely linked with the 
self-regulation connected to feedback processes. 

1.4 The Need for Feedback in Cosmology 

One of the main aims of physical cosmology is to understand in detail 
galaxy formation starting from the initial density fluctuation field with a 
given spectrum (typically a CDM one). Such ab initio computations re- 
quire a tremendous amount of physical processes to be included before it 
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Figure 1.2. Equilibrium solution curve L 
thctical cooling function. See text. 



in the P — T plane for a hypo- 



becomes possible to compare their predictions with experimental data. In 
particular, it is crucial to model the interstellar medium of galaxies, which 
is know to be turbulent, have a multi-phase thermal structure, can undergo 
gravitational instability and form stars. To account for all this complex- 
ity, in addition one should treat correctly all relevant cooling processes, 
radiative and shock heating (let alone magnetic fields!). This has proven 
to be essentially impossible even with present day best supercomputers. 
Hence one has to resort to heuristic models where simplistic prescriptions 
for some of these processes must be adopted. Of course such an approach 
suffers from the fact that a large number of free parameters remains which 
cannot be fixed from first principles. These are essentially contained in 
each of the ingredients used to model galaxy formation, that is, the evolu- 
tion of dark halos, cooling and star formation, chemical enrichment, stellar 
populations. Fortunately, there is a large variety of data against which the 
models can be tested: these data range from the fraction of cooled baryons 
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to cosmic star formation histories, the luminous content of halos, luminos- 
ity functions and faint galaxy counts. The feedback processes enter the 
game as part of such iterative try-and-learn process to which we are bound 
by our ignorance in dealing with complex systems as the galaxies. Still, we 
are far from a full understanding of galaxies in the framework of structure 
formation models. The hope is that feedback can help us to solve some 
"chronic" problems found in cosmological simulations adopting the CDM 
paradigm. Their (partial) list includes: 

(i) Overcooling: the predicted cosmic fraction of cooled baryons is larger 
than observed. Moreover models predict too many faint, low mass 
galaxies; 

(ii) Disk Angular Momentum: the angular momentum loss is too high 
and galactic disk scale lengths are too small; 

(iii) Halo Density Profiles: profiles are centrally too concentrated; 

(iv) Dark Satellites: too many satellites predicted around our Galaxy. 

1.4.1 The Overcooling Problem 

Among the various CDM problems, historically the overcooling has been 
the most prominent and yet unsolved one. In its original formulation it has 
been first spelled out by White & Frenk (1991). Let us assume that, as 
a halo forms, the gas initially relaxes to an isothermal distribution which 
exactly parallels that of the dark matter. The virial theorem then relates 
the gas temperature T-^ir to the circular velocity of the halo Vc, 

kTyir = -^J.mpV^ or T^ir = ■i6V^,.„^/^K, (1.35) 

where iimp is the mean molecular weight of the gas. At each radius in this 
distribution we can then define a cooling time as the ratio of the specific 
energy content to the cooling rate, 

ni{r)A{T) 

where Pg{r) is the gas density profile and rie(r) is the electron density. A(T) 
is the cooling function. The cooling radius is defined as the point where the 
cooling time is equal to the age of the universe, i.e. tcooiifcooi) = t Hubble = 

Hiz)-\ 

Considering the virialized part of the halo to be the region encompass- 
ing a mean overdensity is 200, its radius and mass are defined by 

r„r = O.IH-Hl + z)'/^V,, (1.37) 

M„r -= O.liGHo)-^! + zf^V^. (1.38) 
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Let us distinguish two limiting cases. When TcooI S> r^ir (accretion 
limited case), cooling is so rapid that the infalling gas never comes to hydro- 
static equilibrium. The supply of cold gas for star formation is then limited 
by the infall rate rather than by cooling. The accretion rate is obtained by 
differentiating equation H1.38|l with respect to time and multiplying by the 
fraction of the mass of the universe that remains in gaseous form: 



d ^ 
'di' 



Mace - fgng^0.1{GHor\l + zf^'^V,^ = O.lbfgflgG-Wl (1.39) 



Note that, except for a weak time dependence of the fraction of the initial 
baryon density which remains in gaseous form, fg , this infall rate does not 
depend on redshift. 

In the opposite limit, rcooi 'C ryir (quasi-static case), the accretion 
shock radiates only weakly, a quasi-static atmosphere forms, and the supply 
of cold gas for star formation is regulated by radiative losses near Vcooi ■ A 
simple expression for the inflow rate is 

Mqst = 4:TTpg{rcool)rlool-J^1'cool- (1-40) 

The gas supply rates predicted by equations (|1.39|) and H1.40|) are 
illustrated in Figure [T3I In any particular halo, the rate at which cold gas 
becomes available for star formation is the minimum between Mace and 

Mqst- 

According to the models of Thomas et al. (1986), the bolometric X- 
ray luminosity of the region within the cooling radius of a galactic cooling 
flow is 

Lx = 2.5M,ooiV,^- 

The predictions of this formula are superposed on Figure 11.31 (dashed 
lines). For large circular velocities the mass cooling rates correspond to 
quite substantial X-ray luminosities. 

Integrating the gas supply rates Mcooi over redshift and halo mass 
distribution, we find that for f2(, = 0.04 most of the gas is used before the 
present time. This is unacceptable, since the density contributed by the 
observed stars in galaxies is less than 1% of the critical density. So, the 
star formation results to be too rapid without a regulating process, i.e. 
feedback. 

To solve this puzzle, Larson (1974) proposed that the energy input 
from young stars and supernovae could quench star formation in small 
protogalaxies before more than a small gas fraction has been converted 
into stars. 

Stellar energy input would counteract radiative losses in the cooling 
gas and tend to reduce the supply of gas for further star formation. We can 
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Figure 1.3. Gas infall rate and cooling rates in dark matter lialos as a function of 
circular velocity and redshift. The infall rate {dotted line) is essentially indepen- 
dent of redshift; the cooling rates {solid lines) are given for redshift z = 0, 1, 3, 7 
and 15 {from bottom to top). Dashed lines give present-day X-ray luminosities 
in erg s~^ produced by gas cooling at the given rate in each halo. The pre- 
dicted temperature of this emission is given on the upper abscissa, (a) A cooling 
function for gas of zero metallicity is assumed, (b) A cooling function for metal 
enriched gas. This Figure is taken from White & Frenk (1991). 



imagine that the star formation process is self-regulating in the sense that 
the star formation rate M* takes the value required for heating to balance 
dissipation in the material which docs not form stars. This produces the 
following prescription for the star formation rate: 



M^Vc, z) = e(K)min(Macc, M^st), 
e{V,)^[l + eo{Vo/K)']-'- 



(1.41) 
(1.42) 



For large Vc the available gas turns into stars with high efficiency because 
the energy input is not sufficient to prevent cooling and fragmentation; for 
smaller objects the star formation efficiency e is proportional to V^. The 
assumption of self regulation at small Vc seems plausible because the time 
interval between star formation and energy injection is much shorter than 
either the sound crossing time or the cooling time in the gaseous halos. 
However, other possibilities can be envisaged. For example, Dekel & Silk 
(1986) suggested that supernovae not only would suppress cooling in the 
halo gas but would actually expel it altogether. The conditions leading to 
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such an event are discussed next. 

1.4.2 Dwarf Galaxies: a Feedback Lab 

The observation of dwarf galaxies has shown well-defined correlations be- 
tween their measured properties, and in particular 

(i) Luminosity-Radius L cc R^ with r = 4 

(ii) Luminosity-Metallicity L (x Z^ with z = 5/2 

(iii) Luminosity- Velocity dispersion L oc V^ with v ~ 4 

A simple model can relate the observed scaling parameters r, z, and 
V to each other (Dekel & Silk 1986). Consider a uniform cloud of initial 
mass Mi = Mg + M^ {Mg is the mass of gas driven out of the system, 
and Mi, is the mass in stars) in a sphere of radius Ri which undergoes 
star formation. The metallicity for a constant yield in the instantaneous 
recycling approximation is given by Z = j/ln(l + Mi,/Mg) where y are the 
yields. 

Let us impose simple scaling relations: Mi oc M (gas mass proportional 
to the dark matter mass) , R o:^ Ri , and y ~ T^ (gas loss has no dynamical 
effect). From the observed scaling relation we have L oc R^ (xV^ . Now we 
should write the analogous relations for structure and velocity that hold 
before the removal, i.e. Mi ex R" and Mi ex V'"'^. 

Consider now the case in which the gas is embedded in a dark halo, 
and assume that when it forms stars the mass in gas is proportional to 
the dark matter inside Ri, Mi ex M. If the halo is dominant, the gas loss 
would have no dynamical effect on the stellar system that is left behind, 
so R ~ Ri and V ~ Vi. The relations for structure and velocity that hold 
before the removal give 

L(xi^'^cxy^ M(xR", V^(xM/R, (1.43) 

so that we obtain 

V ^ 2r/(ri - 1). (1.44) 

In the limit M* ^ Mg and L ex M*, from the metallicity relation we 
have 

Z (X :?^ ex — ex Ll/^ (1.45) 

Mg M ^ ^ 

so that 

r/ri = {z-l)/z. (1.46) 
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The final equation conies from tlie energy condition. For tliernial en- 
ergy, in tlie limit of substantial gas loss, Mg ~ Mi, we have L oc MV'^, and 
so 

V = 2z. (1.47) 

For a CDM spectrum P{k) = Ck"^^, we obtain r.^ = 6/(5 + Ug) and 
introducing z = 2r/{r — 1) and Ug = 12/(r + 1) — 5. For r = 4 we find 
Us = —2.6 consistent with the CDM and the scaling relations are 

Loc R'^ oz Z^-^ ex V^-^. (1.48) 

Let us now investigate the critical conditions, in terms of gas density 
n and virial velocity V, for a global supernova-driven gas removal from a 
galaxy while it is forming stars. Here, spherical symmetry and the pres- 
ence of a central point source are assumed. The basic requirement for gas 
removal is that the energy that has been pumped into the gas is enough 
to expel it from the protogalaxy. The energy input in turn depends on the 
supernova rate, on the efficiency of energy transfer into the gas, and on the 
time it takes for the SNRs to overlap and hence affect a substantial fraction 
of the gas. The first is determined by the rate of star formation, the second 
by the evolution of the individual SNRs, and the third by both. When 
all these are expressed as a function of n and V, the critical condition for 
removal takes the form 



E{n,V)>-MgV\ (1.49) 



This relation defines a locus in the n — V diagram shown in Figure 11.51 
within which substantial gas loss is possible. In Figure 11.51 the cooling 
curve, above which the cooling time is less than the free-fall time, confines 
the region where the gas can contract and form stars. The almost vertical 
line Vcrit divides the permissible region for galaxy formation in two; a 
protogalaxy with V > Vcrit would not expel a large fraction of its original 
gas but rather turn most of its original gas into stars to form a 'normal' 
galaxy. A protogalaxy with V < Vcrit can produce a supernova-driven wind 
out of the first burst of star formation, which would drive a substantial 
fraction of the protogalactic gas out, leaving a diffuse dwarf. 

The short-dashed curve marked "Icr" corresponds to density pertur- 
bations SM/M at their equilibrium configuration after a dissipationless 
collapse from a CDM spectrum, normalized to SM/M at a comoving ra- 
dius 8h~^ Mpc. The density n is calculated for a uniformly distributed 
gas in the CDM halos, with a gas-to-total mass ratio x = 0.1- The corre- 
sponding parallel short-dashed curve corresponds to the protogalactic gas 
clouds, after a contraction by a factor x^^ = 10 inside isothermal halos, 
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to densities such that star formation is possible. The vertical dashed ar- 
row marks the largest galaxy that can form out of a typical la peak in 
the initial distribution of density fluctuations. The vast majority of such 
protogalaxies, when they form stars, have V < Vcru, so they would turn 
into dwarfs. The locus where "normal" galaxies are expected to be found 
is the shaded area. It is evident that most of them must originate form 2a 
and 3a peaks in the CDM perturbations. 

So, the theory hence predicts two distinct types of galaxies which oc- 
cupy two distinct loci in the n — V diagram: the "normal" galaxies are 
confined to the region of larger virial velocities and higher densities, and 
they tend to be massive; while the diffuse dwarfs are typically of smaller 
velocities and lower densities, and their mass in star is less than 5x10^ Mq. 
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Figure 1.4. Gas number density vs. virial velocity (or viral temperature), the 
formation of dwarf vs. 'normal' galaxies in CDM halos, and the origin of biased 
galaxy formation. This Figure is taken from Dekel & Silk (1986). 
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This simple model, in spite of its enlightening power, has clear lim- 
itations. The most severe is that is assumes a spherical geometry and a 
single supernova explosion. These assumptions have been released by a 
subsequent study based on a large set of numerical simulations (Mac-Low 
& Ferrara 1999). These authors modelled the effects of repeated supernova 
explosions from starbursts in dwarf galaxies on the interstellar medium of 
these galaxies, taking into account the gravitational potential of their dom- 
inant dark matter halos. They explored supernova rates from one every 
30,000 yr to one every 3 Myr, equivalent to steady mechanical luminosities 
of L = 0.1 — 10 X 10"^* ergs s~^, occurring in dwarf galaxies with gas masses 
Mg = 10^ - 10^ Mq. Surprisingly, Mac-Low & Ferrara (1999) found that 
the mass ejection efficiency is very low for galaxies with mass Mg > 10^ Mq. 
Only galaxies with Mg < IO^Mq have their interstellar gas blown away, 
and then virtually independently of L (see Table n~T|) . On the other hand, 
metals from the supernova ejecta are accelerated to velocities larger than 
the escape speed from the galaxy far more easily than the gas. They found^ 
that for L38 — 1, only about 30% of the metals are retained by a 10^ Mq 
galaxy, and virtually none by smaller galaxies (see Table IL^ . 



Visible Mass 


Mechanical Luminosity [10^* erg s ^] 


[Afo] 


0.1 1.0 10 



10« 0.18 1.0 1.0 

lO'' 3.5 X 10"^ 8.4 X 10"^ 4.8 x 10"^ 



-5 



10^ 1.1 X 10"" 3.4 X 10"* 1.3 X 10"^ 

10^ 0.0 7.6 X lO"'* 1.9 X 10 



Table 1.1. Mass ejection efficiency 
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10^ 
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1.0 


10* 


0.8 


1.0 


1.0 


10^ 


0.0 


0.69 


0.97 



Table 1.2. Metal ejection efficiency 



^ Throughout the text we use the standard notation Yx = Y/IQ-^ 
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Figure 1.5. Distribution of metal-enriclied stellar outflow material and SN ejecta 
from the starburst energy source at time 200 Myr. Note that in most cases, no 
enriched gas remains in the disks of the galaxies. This Figure is taken from 
Mac-Low & Ferrara (1999). 



1.4.3 Blowout, Blowaway and Galactic Fountains 

The results of the MacLow & Ferrara study served to clearly classify the 
various events induced by starburst in galaxies. We can distinguish between 
blowout and blowaway processes, depending on the fraction of the parent 
galaxy mass involved in the mass loss phenomena. Whereas in the blowout 
process the fraction of mass involved is the one contained in cavities created 
by the supernova or superbubble explosions, the blowaway is much more 
destructive, resulting in the complete expulsion of the gas content of the 
galaxy. The two processes lie in different regions of the (1 + ^) — Mg^i 
plane shown in figure H"^ where (j> = Mh/Mg is the dark-to- visible mass 
ratio and Mgj = Mg/10^ Mq (Ferrara & Tolstoy 2000). Galaxies with 
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gas mass content larger than 10^ Mq do not suffer mass losses, due to 
their large gravitational well. Of course, this does not rule out the possible 
presence of outflows with velocities below the escape velocity (fountain) 
in which material is temporarily stored in the halo and then returns to 
the main body of the galaxy. For galaxies with gas mass lower than this 
value, outflows cannot be prevented. If the mass is reduced further, and 
for ^ < 20, a blowavay, and therefore a complete stripping of the galactic 
gas, should occur. To exemplify, the expected value of (j) as function of Mg 
empirically derived by Persic et al. (1996) is also plotted, which should 
give an idea of a likely location of the various galaxies in the (1 + ^) ~ Mgj 
plane. 

14.3.1 Blowout 

The evolution of a point explosion in an exponentially stratified medium 
can be obtained from dimensional analysis. Suppose that the gas density 
distribution is horizontally homogeneous and that p{z) — pQexp{—z/H). 
The velocity of the shock wave is w ~ (P/p)^^^, where the pressure P is 
roughly equal to E/z^, and E is the total energy if the explosion. Then it 
follows that 

v{z) ~ -B'/'po '^' exp(z/2iJ)z-3/2. ^^ ^q^ 

This curve has a minimum at z = 3H and this value defines the height at 
which the shock wave, initially decelerating, is accelerated to infinity and 
a blowout takes place. Therefore, 3H can be used as the fiducial height 
where the velocity Vb = v(3H) is evaluated. 

There are three different possible fates for SN-shocked gas, depending 
on the value of Vb. If Vb < Cs^eff, where Cg is the sound speed in the ISM, 
then the explosion will be confined in the disk and no mass-loss will occur; 
for Cs,eff < Vb < Ve (ve IS the escape velocity) the supershell will breakout 
of the disk into the halo, but the flow will remain bound to the galaxy; 
finally, Ve < Ufc will lead to a true mass-loss from the galaxy. 

1.4-3.2 Blowaway 

The requirement for blowaway is that the momentum of the shell is larger 
than the momentum necessary to accelerate the gas outside the shell at 
velocity larger than the escape velocity, MgVc < MoVe- Defining the disk 
axis ratio as e = R/H > 1 , the blowaway condition can be rewritten as 

^>(e-a)2a-V/2, (1.51) 



Ve 



where a — 2/3. The above equation is graphically displayed in Figure [TTTI 
Flatter galaxies (large e values) preferentially undergo blowout, whereas 
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Figure 1.6. Regions in the (1 + 0) — Mg,7 plane in which different dynamical 
phenomena may occur. Also shown is the locus point describing the Persic et al. 
(PSS) relation (dotted line). This Figure is taken from Ferrara & Tolstoy (2000). 



rounder ones are more likely to be blown-away; as Vb/ve is increased the 
critical value of e increases accordingly. Unless the galaxy is perfectly 
spherical, blowaway is always preceded by blowout; between the two events 
the aspect of the galaxy may look extremely perturbed, with one or more 
huge cavities left after blowout. 
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Figure 1.7. Conditions for blowaway, blowout and confinement as a function of 
tiie major-to-minor axis ratio e of dwarf galaxies; e — 1 corresponds to spherical 
bodies. This Figure is taken from Ferrara & Tolstoy (2000). 



1.4.4 Further Model Improvements 

All models discussed so far make the so-called SEX BOMB assumption, 
i.e. a Spherically Expanding blastwave. This is a good assumption as long 
as a single burst regions exists whose size is small compared to the size of 
the system and it is located at the center of the galaxy. This however is a 
rather idealized situation. 

A more detailed study about the possibility of mass loss due to dis- 
tributed SN explosions at high redshift has been carried on by Mori, Ferrara 
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& Madau (2002). These authors presented results from three-dimensional 
numerical simulations of the dynamics of SN-driven bubbles as they prop- 
agate through and escape the grasp of subgalactic halos with masses M = 
10^ h~^ Mq at redshift z = 9. Halos in this mass range are characterized 
by very short dynamical timescales (and even shorter gas cooling times) 
and may therefore form stars in a rapid but intense burst before SN 'feed- 
back' quenches further star formation. This hydrodynamic simulations use 
a nested grid method to follow the evolution of explosive multi-SN events 
operating on the characteristic timescale of a few x 10^ yr, the lifetime of 
massive stars. The results confirm that, if the star formation efficiency 
of subgalactic halos is « 10%, a significant fraction of the halo gas will 
be lifted out of the potential well ('blow-away'), shock the intergalactic 
medium, and pollute it with metal-enriched material, a scenario recently 
advocated by Madau, Ferrara, & Rees (2001). Depending on the stellar 
distribution, Mori et al. (2002) found that less than 30% of the available 
SN energy gets converted into kinetic energy of the blown away material, 
the remainder being radiated away. 

However, it appears that realistic models lead to the conclusion that 
mechanical feedback is less efficient than expected from SEX BOMB sim- 
ple schemes. The reason is that off-nuclear SN explosions drive inward- 
propagating shocks that tend to collect and pile up cold gas in the central 
regions of the host halo. Low-mass galaxies at early epochs then may 
survive multiple SN events and continue forming stars. 

Figure 11.81 and 11.91 show the fraction of the initial halo baryonic mass 
contained inside (1, 0.5, 0.1) of the virial radius r^ir as a function of time 
for two different runs: an extended stellar distribution (case 1) and a more 
concentrated one (case 2). The differences are striking: in case 1, the 
amount of gas at the center is constantly increasing, finally collecting inside 
0.1r„ir about 30% of the total initial mass. On the contrary, in case 2, the 
central regions remain practically devoided of gas until 60 Myr, when the 
accretion process starts. The final result is a small core containing a fraction 
of only 5% of the initial mass. In the former case, 50% of the halo mass 
is ejected together with the shell, whereas in case 2 this fraction is ^ 85%, 
i.e., the blow-away is nearly complete. 

As a first conclusion, it seems that quenching star formation in galaxies 
by ejecting large fraction of their gas seems very difficult and hence unvi- 
able as a feedback scheme (although this might be possible in very small 
galaxies) . Heavy elements can instead escape much more easily as they are 
carried away by the mass-unloaded, hot, SN-shocked gas; energy is carried 
away efficiently as well. This conclusion raises the issue if star formation 
must be rather governed by more gentle and self-regulating processes, i.e. 
a 'true' feedback. 
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Figure 1.8. Snapshots of the logarithmic number density of the gas at five 
different elapsed times. The three panels in each row show the spatial den- 
sity distribution in the x — y plane on the nested grids. The density range is 
-5 < log(n/cm"^) < 1. This Figure is taken from Mori et al. (2002) 
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Figure 1.9. The evolution of the gas mass inside the gravitational potential well 
of the CDM halo for an extended stellar distribution (case 1) as a function of 
time. Curves correspond to the gas mass inside the virial radius rvir {solid line), 
O.Srvir {dashed line), and O.lrvir {dotted line). This Figure is taken from Mori et 
al. (2002). 

1.5 The Gentle Feedback 



In numerical simulations an astonishingly large number of attempts to im- 
plement feedback scheme have been tried. In spite of this wealth, essen- 
tially all of them can be classified in one of the two following categories: 
(i) 'Thermal' feedback schemes in which the energy released by a super- 
nova is assumed to simply heat the ISM. The dynamics of the gas is only 
slightly affected: this is because the dense gas in star-forming region is able 
to radiate away this heat input very quickly. Moreover it would imply a 
huge amount of hot gas, that is hardly seen in local universe, (ii) 'Kinetic' 
feedback schemes, in which the explosion energy is assigned to fluid parcels 
as pure kinetic energy. This algorithm performs slightly better and it is 
probably more physical. It would be impossible to review in detail here 
even only a few of these schemes. We have thus chosen to describe a min- 
imum feedback model that tries to catch the most important feature that 
is missing in the previous ones, that is the multiphase structure and the 
matter exchange among gas phases characterizing real galaxies. 
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Figure 1.10. The same of Figure [T!^ but for a more concentrated stellar distri- 
bution (case 2). This Figure is taken from Mori et al. (2002). 



1.5.1 Feedback as ISM Sterilization 

Dynamical models of the interstellar medium (ISM) describe this complex 
system in terms of a dense, cold neutral component (CNM) co-spatial with a 
warm neutral intercloud medium (WNM) by which it is pressure confined. 
Such a multiphase medium offers an efficient self-regulatory mechanism 
for star formation. Schematically, as more stars are formed, the extra 
heating caused by the enhancement of their UV field tends to transfer mass 
from the CNM into the WNM, where the conditions are highly unfavorable 
to star formation, i.e. a sterile phase. This acts to decrease the star 
formation rate and slowly brings back gas from the WNM reservoir in the 
cold phase where stars can start to form again. Under most conditions 
this cycle is stable and tends to regulate the amount of gas turned in stars 
quite effectively. In this scenario starburst can be only produced either by 
external triggers (mergings and interactions with other galaxies) and/or 
during the very first phase of galaxy formation, when the gentle feedback 
has not yet had the time to control the system. The physical basis of the 
gentle feedback is the thermal instability already discussed early on. Field 
(1965) discussed the criteria for stability of gas in thermal equilibrium (i.e. 
L{n,T) = nT{n,T) - n^k{n,T) = 0, see also sect. [Oil. In a plot of the 
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thermal pressure P/k versus the hydrogen density n, the region of thermal 
stability occur for (i(logP)/(i(logn) > 0. Thus, Figure IT . 1 II shows that a 
stable two-phase medium, can be maintained in pressure equilibrium for 
Pmin < P < Pmax- For pressurc less than Pmin only the warm phase 
is possible, while at pressures greater than Pmax only the cold phase is 
possible (Wolfire et al. 1995). Moreover, in Figure IT Til it is shown, that 
when the metallicity of the gas Z, assumed to be equal to the dust-to-gas 
ratio D/G, is decreased, the pressure range in which a multiphase medium 
can exist becomes wider; in addition, the mean equilibrium pressure and 
density of the cold gas increase (Ricotti et al. 1997). 

The first advantage of such feedback scenario is that there is no need 
for gas loss as required by canonical feedback schemes: typically, the loss 
rate necessary to regulate the star formation activity is difficult to sustain 
for galaxies with mass of the order of the Milky Way for interestingly long 
periods. In addition, the production of large amounts of hot, X-ray emitting 
gas is expected. Such gas is rarely seen in/around observed galaxies. 

An additional advantage of these models is the self-regulatory behavior 
of the process, that we have seen to be typical of feedback. In fact, there is 
no need of infall or outflow of gas from the galaxy (close box). As cold gas 
is available, the galaxy has a star formation burst, leading to the increase of 
the metallicity and of the dust content. Moreover, extra heating is provided 
by the SN explosions and from UV background, so that some gas mass is 
transferred from the CNM to the WIM. The star formation is quenched 
since the gas is in the warm, sterile phase, and the heating drops down. 
After a cooling time, the gas becomes again available for star formation. 

The disadvantage of these models is that they are physically complex, 
since the involve the interplay between dynamics and thermodynamics. 
Hence, a number of physical processes await inclusion in current cosmolog- 
ical simulations. 

1.5.2 A Porosity- regulated Feedback 

A similar feedback scheme has been recently proposed by Silk (2003). In 
such a feedback prescription, the filling factor of the hot, SN-shocked gas 
plays a central role. The filling factor fhot of hot gas can be expressed in 
terms of the porosity Q as 

hot^l-e'Q, (1.52) 

where Q is regulated by the pressure confinement of the supernova bubbles. 
As fhot — *■ 1, the star formation rate reaches a saturation level and for / = 1 
a 'buoyant' wind is produced with outflow rate approximately equal to the 
star formation rate. 

This feedback has the advantage to be based on the physical treatment 
of the supernova explosion, that current numerical simulations fail to follow 
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Figure 1.11. Thermal pressure P/k vs. hydrogen density n; the curves refer to 
different values of the dust-to-gas ratio, D/G, and metallicity Z, with D/G = Z. 
The gas is thermally stable for (d log P/d log n) > (i.e. positive slope of the 
curves). Unless Z — 0, a stable two-phase medium is supported. This Figure is 
taken from Ricotti et al. (1997). 

in detail due to the lack of resolution. Moreover the expression of the 
porosity Q is identical for all the masses: the outflow rate is of the order 
of the star formation rate. Only the star formation efficiency depends on 
potential well depth {ecr oc V'j^'^). The wind efficiency depends only weakly 
on 14. In summary, outflows can occur even for massive galaxies since star 
formation efficiency is larger in these systems. 

However, there are some problems. The porosity-driven feedback is 
not a true multiphase model and has to be generalized to include the time 
evolution of Q, metallicity, UV background, etc. Moreover, it again implies 
a large amount of hot gas that is not detected in galaxies. 



1.5.3 Advanced Multiphase/Feedback Schemes 

The main problem affecting numerical simulations of feedback processes is 
that, owing to the limited numerical resolution, there is a large disparity 
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between the niinimuni simulation timescale and the true physical timescale 
associated with supernova feedback (Thacker & Couchman 2001). Phys- 
ically, following a supernova event, the gas cooling time very rapidly in- 
creases to O(IO^) yr over the time it takes the shock front to propagate. 
Thus, it is the timescale - t^ - for the cooling time to increase markedly 
that is of critical importance. In the SPH simulations, 'stars' form and 'su- 
pernova feedback' occurs in gas cores that are typically less than /imm/5 in 
diameter {hmin being the effective spatial resolution). Any rearrangement 
of the particles in such a small region leads to a density estimate varying 
by at most ~ 7%. Since the cooling time in the model is dependent on 
the local SPH density, any energy deposited into a feedback region only 
reduces the local cooling rate by increasing the temperature; the SPH den- 
sity cannot respond on a timescale that in real supernova events would 
very rapidly increase the cooling time. A successful feedback algorithm in 
an SPH model must thus overcome the fact that psPH does not change on 
the same timescale as i^ and consequently that cooling times for dense gas 
cores remain short. 

Standard implementations of SPH overestimate the density for parti- 
cles that fall near the boundary of a higher density phase (Pearce et al. 
1999). The usual assumption of small density gradients across the smooth- 
ing kernel breaks down in this regime, and nearby clusters of high density 
particles cause an upward bias in the standard SPH estimate: 

N 

< Pi>^^mjW^j, (1.53) 

j 

where N is the number of neighbors j of particle i and Wij is a symmetric 
smoothing kernel. In order to avoid this bias, which leads to artificially 
high cooling rates and to spatial exclusion effects, the neighbor search and 
the density evaluation of equation H1.53|l has to be modified in a way which 
leaves the numerical scheme as simple as possible. It is important to recog- 
nize that the local density is the gas property responsible for phase segre- 
gation (since it determines the local cooling rate). An appropriate density 
estimator for a particle of any given phase should use only local material 
which is also part of that same phase. 

Marri & White (2003) proposed a new feedback scheme that takes into 
account these considerations. This works as follows. As a supernova ex- 
plodes, the energy is distributed in a fraction €c (e^,) to the cold (warm) gas. 
The fundamental hypothesis is that at large scale (of kpc order, say) the net 
effect of all the complex microscopic processes is well described by an en- 
ergy input shared by the macroscopic phases in given proportions. Values 
for Ec and e™ could be fixed from a complete theory of the ISM, describing 
all the relevant processes, or through numerical simulations. Feedback to 
the hot phase is implemented by adding thermal energy to the ten near- 
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est hot neighbors. Feedback to the cold phase is instead accumulated in 
a reservoir within the star-forming particle itself. This continues until the 
accumulated energy is sufficient to heat the gas component of the particle 
above T ^ SOT* ~ 10^ K. This is far enough above T^ for the promoted 
particle to be considered 'hot' in its subsequent evolution. At the same 
time any 'hidden' stellar content is dumped to the cold neighbors. 

A simple test of these new scheme (the initial condition are the ro- 
tating, centrally concentrated sphere described in Navarro & White (1993) 
which consists of 90% by mass in dark matter) shows that hot particle do 
survive near the central disk in the multiphase SPH case. They have an 
almost spherically symmetric distribution with density peaked at the cen- 
ter f Figure ri.l2|l . The cold disk rotates within this ambient hot medium. 
In contrast, in the standard SPH model hot particle are excluded from the 
vicinity of the disk so that the hot phase actually has a density minimum 
at the center of the galaxy. This is seen most clearly in the gas density 
profile of Figure 11.131 Such behavior is clearly unphysical. 



1.6 Feedbacks in the Early Universe 

At z = 1000 the Universe has cooled down to 3000 K and the hydrogen 
becomes neutral (recombination). Then, at z < 20 the first stars (Popula- 
tion or Pop III stars) form and these gradually photo-ionize the hydrogen 
in the IGM (reionization) . These epoch, witnessing the return of light in 
the universe after the Big Bang, is usually dubbed as the Cosmic Dawn. 
At z < 6, galaxies form most of their stars and grow by merging. Finally 
at z < 1, the massive clusters are assembled. 

We can distinguish three main feedback types that are fundamentally 
shaping the universe at cosmic dawn: 

(i) Stellar Feedback: the process of star formation produces the con- 
ditions (destruction of cold gas) such that astration (or even galaxy 
formation) comes to a halt or is temporarily blocked; 

(ii) Chemical Feedback: massive stars by polluting the ambient gas 
with heavies, shift the prevailing formation mode toward low mass 
stars (self-killing population) 

(iii) Radiative Feedback: the process of star formation erases the con- 
ditions (molecular formation) necessary to form stars via UV radia- 
tion and ionizing photons. Such radiative input heats and ionizes the 
gas decreasing its ability to produce sources (stars) to further sustain 
reionization. 

In the following we will describe in detail the three feedback types. 
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Figure 1.12. Face-on projections at i ~ 1.2 Gyr of 8000 gas particle, cooling 
only simulations of a collapsing, rotating sphere. The standard SPH simulation 
is in the upper row with the MSPH simulation below it. The left-hand plots 
show "hot" particles (T > 10^ K) while the right hand ones show "cold" particles 
(r < lO^K and uh > 0.1 cm"^). This Figure is taken from Marri & White 
(2003). 



1.6.1 Stellar Feedback 

1.6.1.1 Oxygen in the Low Density IGM 

One of the most obvious signatures of stellar feedback is the metal en- 
richment of the intergalactic medium. If powerful winds are driven by 
supernova explosions, one would expect to see widespread traces of heavy 
elements away from their production sites, i.e. galaxies. Schaye et al. 
(2000) have reported the detection of OVI in the low-density IGM at high 
redshift. They perform a pixel- by-pixel search for OVI absorption in eight 
high quality quasar spectra spanning the redshift range z = 2.0 — 4.5. At 
2 < z < 3 they detect OVI in the form of a positive correlation between 
the HI Lya optical depth and the optical depth in the corresponding OVI 
pixel, down to thi ^0.1 (underdense regions), that means that metals 
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Figure 1.13. Particle densities as a function of galactocentric distance after 
~ 1.2 Gyr of evolution in two 8000 gas particle, cooling only simulations of the 
collapse of a rotating sphere. The upper plot is for an SPH model and the lower 
for an MSPH model. Gas particles with T > lO^K are plotted with crosses. This 
Figure is taken from Marri & White (2003). 



are far from galaxies (Figure ri.l4|l . Moreover, the observed narrow widths 
of metal absorption lines (CIV, SilV) lines lines imply low temperatures 
TiGM '^ few X lO-* K. 

A natural hypothesis would be that the Lya forest has been enriched 
by metals ejected by Lyman Break Galaxies at moderate redshift. The 
density of these objects is riLBG = 0.013 h^ Mpc'^. A fiUing factor of 
~ 1% is obtained for a shock radius Rg — 140 h~^ kpc, that corresponds 
at z — 3 {h — 0.5) to a shock velocity Vg — 600 km s^^. In this case, 
we expect a post shock gas temperature larger than 2 x 10^ K, that is 
around hundred times what we observed. So the metal pollution must have 
occurred earlier than redshift 3, resulting in a more uniform distribution 
and thus enriching vast regions of the intergalactic space. This allows the 
Lya forest to be hydrodynamically 'cold' at low redshift, as intergalactic 
baryons have enough time to relax again under the influence of dark matter 
gravity only (Scannapieco, Ferrara & Madau 2002). 
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In Figure 11.151 is shown the thermal history of the IGM as a func- 
tion of redshift as computed by Madau, Ferrara & Rees (2001). The gas 
is ahowed to interact with the CMB through Compton coohng and either 
with a time-dependent quasar-ionizing background as computed by Haardt 
& Madau (1996) or with a time-dependent metagalactic flux of intensity 
10~^^ ergcm^^ s^^ Hz~^ sr^^ at 1 Ryd and power-law spectrum with en- 
ergy slope a — I. The temperature of the medium at 2: = 9 has been either 
computed self-consistently from photoheating or fixed to be in the range 
lO'*'^ — 10^ K, as expected in SN-driven bubbles with significant filling fac- 
tors. The various curves show that the temperature of the IGM at z = 3 — 4 
will retain little memory of an early era of pregalactic outflows. 

The large increase of the IGM temperature at high redshift connected 
with such an era of pregalactic outflows, causes a larger Jeans mass, thereby 
preventing gas from accreting efficiently into small dark matter halos (Ben- 
son & Madau 2003). For typical preheating energies, the IGM is driven to 
temperatures just below the virial temperature of halos hosting L^, galaxies. 
Thus we may expect preheating to have a strong effect on the galaxy lumi- 
nosity function at 2; = 0. Moderate preheating scenarios, with Tiqm ^ 10^ 
K at z ~ 10, are able to flatten the faint-end slope of the luminosity func- 
tion, producing excellent agreement with observations, without the need 
for any local strong feedback within galaxies (Figure [1.1611 . 

1.6.1.2 Gas Stripping from Neighbor Galaxy Shocks 

The formation of a galaxy can be inhibited also by the outflows from neigh- 
boring dwarfs as the result of two different mechanisms (Scannapicco et al. 
2000). In the "mechanical evaporation" scenario, the gas associated with 
an overdense region is heated by a shock above its virial temperature. The 
thermal pressure of the gas then overcomes the dark matter potential and 
the gas expands out of the halo, preventing galaxy formation. In this case, 
the cooling time of the collapsing cloud must be shorter than its sound 
crossing time, otherwise the gas will cool before it expands out of the grav- 
itational well and will continue to collapse. 

Alternatively, the gas may be stripped from a collapsing perturbation 
by a shock from a nearby source. In this case, the momentum of the shock 
is sufficient to carry with it the gas associated with the neighbor, thus 
emptying the halo of its baryons and preventing a galaxy from forming. 

In principle, outflows can suppress the formation of nearby galaxies 
both by shock heating/evaporation and by stripping of the baryonic matter 
from collapsing dark matter halos; in practice, the short cooling times 
for most dwarf-scale collapsing objects suggest that the baryonic stripping 
scenario is almost always dominant f Figure fl. 171) . This mechanism has the 
largest impact in forming dwarfs in the < lG^Af0 range which is sufficiently 
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Figure 1.14. IGM Lya-OVI optical depth correlation from the pixel analysis 
of the spectra of 8 QSOs. Vertical error bars are la errors, horizontal error bars 
are smaller than the symbols and are not shown. For z < 3 Tovi.app and thi are 
clearly correlated, down to optical depths as low as thi ~ 10^^. A correlation 
between rovi.app and thi implies that OVI absorption has been detected in the 
Lya forest (Schaye et al. 2000). 



large to resist photoevaporation by UV radiation, but too small to avoid 
being swept up by nearby dwarf outflows. 



1.6.2 Chemical Feedback 

Recent studies (Bronim, Coppi & Larson 1999, 2002; Abel, Bryan & Nor- 
man 2000; Schneder et al. 2002; Omukai & Inutsuka 2002; Omukai & Palla 
2003) have shown that in absence of heavy elements the formation of stars 
with masses 100 times that of the Sun would have been strongly favored, 
and that low-mass stars could not have formed before a minimum level of 
metal enrichment had been reached. 

Bromm et al. (2001) have studied the effect of the metallicity on the 
evolution of the gas in a collapsing dark matter mini-halo, simulating an 
isolated 3<t peak of mass 2 x 10^ M© that collapse at 2; ~ 30, using smoothed 
particle hydrodynamics. The gas has a supposed level of pre-enrichment 
of either Z = 10^^ Zq of lO^'' Zq. The H2 is assumed to be radiatively 
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Figure 1.15. Thermal history of intergalactic gas at the mean density in an 
Einstein-de Sitter universe with Qth^ ~ 0.019 and h = 0.5. Short dash-dotted 
line: temperature evolution when the only heating source is a constant ultravi- 
olet (CUV) background of intensity 10~^^ ergcm"^ s"'^ Hz"'^ sr~^ at 1 Ryd and 
power-law spectrum with energy slope a = 1. Long dash-dotted line: same 
for the time-dependent quasar ionizing background as computed by Haardt & 
Madau (1996; HM). Short dashed line: heating due to a CUV background but 
with an initial temperature of 4 x 10'' K at z = 9 as expected from an early era 
of pregalactic outflows. Long dashed line: same but for a HM background. Solid 
Itne: heating due to a CUV background but with an initial temperature of lO"" K 
at 2 = 9. Dotted line: same but for a HM background. This Figure is taken from 
Madau, Ferrara & Rees (2001) 



Feedbacks in the Early Universe 



35 



10° 



10 



10- 



o 






-& 



10- 



10- 



10- 




No preheating 



■/ ^ = 9, E ^ =0.1keV 

preheat ' preheat 

Veheal"12, Ep^^^^^^- . 1 keV 

■/ , ,= 12. E , , = 0.3koV 

preheat preheat 



Norberg et al. (2002) 



-15 



-20 



M^ -51ogh 



Figure 1.16. B-band luminosity functions of galaxies at z — 0, as predicted by 
the semi-analytic model of Benson et al. (2002). The observational determination 
of Norberg et al. (2002) is shown as circles. This Figure is taken from Benson & 
Madau (2003). 



destroyed by the presence of a soft UV background. Moreover Bronim et 
al. (2001) do not consider the presence of molecules or dust. 

The evolution proceeds very differently for the two cases. The gas in 
the lower metallicity case fails to undergo continued collapse and fragmen- 
tation, remaining in the post-virialization, pressure-supported state in a 
roughly spherical configuration. The final result is likely to be the forma- 
tion of a single very massive (~ 10"^ — 10'' Mq) object (Figure ri.l8|l . On 
the other hand fFigure [1.19(l . the gas in the higher metallicity case can cool 



Feedbacks in the Early Universe 



36 



10 



^ 







9 



8 




Pho toe vapor at ion 



7 



Mech Evap_ 




I Photoevaporation" 







5 



10 



15 



Figure 1.17. Relevant mass scales for suppression of dwarf galaxy formation. 
The upper lines are the masses below which halos will be heated beyond their 
virial temperatures, although cooling prevents mechanical evaporation from oc- 
curring for halos with masses above the cross hatched regions. The second high- 
est set of lines, bounding the lightly shaded regions, show the masses below 
which baryonic stripping is effective. Finally the heavily shaded regions show 
objects that are susceptible to photoevaporation. The upper panel is a flat CDM 
model and the lower panel is a flat ACDM model with fio ~ 0.3. In all cases 
(eNh) = 5000r2|^\ the overdensity 5 = p/po = 2.0, fit = 0.05, h = .65. Note 
that photoevaporation affects a larger mass range than mechanical evaporation 
in the ACDM cosmology. This Figure is taken from Scannapieco et al. (2000) 
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efficiently and collapse into a disk-like configuration. The disk material is 
gravitational unstable and fragments into a large number of high density 
clumps, that are the seeds of low mass stars. So, they conclude that there 
exists a critical metallicity Zcrit at which the transition from the formation 
of high mass objects to a low-mass star formation mode, as we observe at 
present, occurs. In the following we refer to Pop III stars in the mass range 
100-600 Mq forming out of the collapse oi Z < Zcrit gas clouds. 
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Figure 1.18. Simulation of a collapsing dark matter halo with mean metallicity 
oi Z — 10^'* Zq. Morphology at 2; = 22.7. Top row; the remaining gas in the 
diffuse phase. Bottom row: distribution of clumps. Dot sizes are proportional to 
the mass of the clumps. Left panels: face-on view. Right panels: edge-on view. 
The box size is 30 pc. It can be seen that no fragmentation occurs. This Figure 
is taken from Bromm et al. (2001). 



More recently, Schneider et al. (2002) have investigated the problem 
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Figure 1.19. Simulation of a collapsing dark matter halo with mean metallicity 
of Z = 10"^ Zq. Morphology at ^ = 28. The convention of Figure [TTgl is 
adopted for the row and columns. It can be seen that the gas has undergone 
vigorous fragmentation. This Figure is taken from Bromm et al. (2001) 



more in detail including the cooling due to H2, molecular, and dust, using 
the model of Omukai (2001). The gas within the dark matter halo is given 
an initial temperature of 100 K, and the subsequent thermal and chemical 
evolution of the gravitationally collapsing cloud is followed numerically un- 
til a central protostellar core forms. The results for different metallicities 
are summarized in Figure 11.201 showing the temperature and adiabatic in- 
dex evolution as a function of the hydrogen number density of protostellar 
clouds for different metalhcity {Z = 0, 10"*^, lO"", 10-^ 1 Zq). The neces- 
sary conditions to stop fragmentation and start gravitational contraction 
within each clump arc that cooling becomes inefficient, and the Jeans mass 
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of the fragments does not decrease any further, thus favoring fragmentation 
into sub-clumps. This condition depends somewhat on the geometry of the 
fragments, and translates into 7 > 4/3 (7 > 1) for spherical (filamentary) 
clumps, where 7 is the adiabatic index defined as T ex nP'^^. 

For a metal-free gas, the only efficient coolant is molecular hydro- 
gen. Cooling due to molecular line emission becomes inefficient at densities 
above n > 10'^ cm^"', and fragmentation stops when the minimum fragment 
mass is of order 10'^— 10^ Mq. Clouds with mean metallicity Z = 10^® Zq 
follow the same evolution as that of the gas with primordial composition in 
the (n, T) plane. For metallicities larger than 10~^ Zq the fragmentation 
proceeds further until the density is ~ 10^^ cm""^ and the corresponding 
Jeans mass is of the order 10~^ Mq. Schneider et al. (2002) concluded 
that the critical metallicity locates in the range 10~^ — 10"'* Zq. 

According to the scenario proposed above, the first stars that form 
out of gas of primordial composition tend to be very massive, with masses 
~ 10^ — 10"^ Mq. It is only when metals change the composition of the gas 
that further fragmentation occurs, producing stars with significantly lower 
masses. A series of numerical studies (Heger & Woosley 2001; Fryer et al. 
2001; Umeda & Nomoto 2002) have investigated the nucleosynthesis and 
final state of metal-free massive stars. Heger & Woosley (2001) dehneate 
three mass ranges characterized by distinct evolutionary paths: 

(i) Mi, > 260 Mq. - The nuclear energy released from the collapse of 
stars in this mass range is insufficient to reverse the implosion. The 
final result is a very massive black hole (VMBH) locking up all heavy 
elements produced. 

(ii) 140 Mq < M* < 260 M©. - The mass regime of the pair-unstable 
supernovae (SN^-y). Precollapse winds and pulsations result in little 
mass loss; the star implodes to a maximum temperature that depends 
on its mass and then explodes, leaving no remnant. The explosion 
expels metals into the surrounding ambient ISM. 

(ui) 30 MQ<Mi,< 140 Mq. - Black hole formation is the most hkely out- 
come, because either a successful outgoing shock fails to occur or the 
shock is so weak that the fallback converts the neutron star remnant 
into a black hole (Fryer 1999) 

Stars that form in the mass ranges 1 and 3 above fail to eject most of 
their heavy elements. If the first stars have masses in excess of 260 Mq, 
they invariably end their lives as VMBHs and do not release any of their 
synthesized heavy elements. However, as long as the gas remains metal 
free, the subsequent generations of stars will continue to be top-heavy. 
This 'star formation conundrum' can be solved only if a fraction /^^ of the 
first generation of massive stars are in the SN-y^ range and enrich the gas 
with heavy elements up to a mean metallicity oi Z > 10~^^^ Zq. 

Schneider et al. (2002) computed the transition redshift Zf at which 
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the mean metallicity is 10"''' Zq for various values of /-y^. The results 
are plotted in Figure 11.211 along with the corresponding critical density 
^VMBH contributed by the VMBHs formed. In order to set some limits on 
fry-y, we have to compare the predicted critical density of VMBH remnants 
to present observational data. This depends on the assumptions about the 
fate of these VMBHs at late times. 

Under the hypothesis that all VMBHs have, during the course of 
galaxy mergers, been used to build up the supermassive black holes (SMBHs) 
detected today so that ^vmbh = ^smbh, from the top panel of Figure 
ll.21l we can derive /^^ ~ 0.06 and Zf ~ 18.5. 

At the other extreme, wherein the assembled SMBHs in galactic cen- 
ters have formed primarily via accretion and are unrelated to VMBHs, two 
possibilities can be distinguish: (1) VMBHs would still be in the process of 
spiraling into the center of galaxies because of dynamical friction but are 
unlikely to have reached the center within a Hubble time because of the 
long dynamical friction timescale (Madau & Rees 2001), or, (2) VMBHs 
contribute the entire baryonic dark matter in galactic halos. In case (1) 
some fraction of these VMBHs might appear as off-center accreting sources 
that show up in hard X-ray wave band. In fact, ROSAT and Chandra have 
detected such objects. Using the contribution to Vl from X-ray-bright, off- 
center i?OS'^T sources as a limit for the density of VMBHs, from the lower 
dashed line in the bottom panel of Figure 11.211 f.y^ has to be nearly unity 
and Zf > 22.1. In case (2), assuming that the baryonic dark matter in 
galaxy halos is entirely contributed by VMBH (upper dashed line in the 
bottom panel of Figure IT. 2 1|) . /^^ ~ 3.15 x 10~^ and Zf > 5.4. Therefore 
the case in which SMBHs are unrelated to VMBHs, only weak limits for 
f^Yj and Zf can be obtained. Although the actual data do not allow us at 
present to strongly constrain these two quantities, they provide interesting 
bounds on the proposed scenario. 

We have seen that there exists a critical metallicity Zcr = 10^^^ Zq 
(Schneider et al. 2002; 2003) at which we expect the transition between 
a top-heavy and a normal Salpeter-like IMF. One point that might have 
important observational consequences is the fact that cosmic metal enrich- 
ment has proceeded very inhomogeneously (Scannapieco, Ferrara & Madau 
2002; Furlanetto & Loeb 2003; Marri et al. 2003), with regions close to 
star formation sites rapidly becoming metal-polluted and overshooting Zcr, 
and others remaining essentially metal-free. Thus, the two modes of star 
formation, Pop HI and normal, must have been active at the same time 
and possibly down to relatively low redshifts, opening up the possibility of 
detecting Pop III stars. 

Scannapieco, Schneider & Ferrara (2003) have studied using an ana- 
lytic model of inhomogcneous structure formation, the evolution of Pop III 
objects as a function of the star formation efficiency, IMF, and efficiency 
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Figure 1.20. Top: Evolution of the temperature as a function of tlie liydrogen 
number density of protostellar clouds with the same initial gas temperature but 
varying the metallicities Z = (0, 10~®, 10~*, 10"'^, 1) Zq [Z increasing from top 
to bottom curves). The dashed lines correspond to the constant Jeans mass for 
spherical clumps; open circles indicates the points where fragmentation stops; 
filled circles mark the formation of hydrostatic cores. Bottom: The adiabatic 
index 7 as a function of the hydrogen number density for the curves shown in 
the top panel. Dotted (dashed) lines correspond to 7 = 1 (7 = 4/3); open and 
filled circles as above. This Figure is taken from Schneider et al. (2002). 



of outflow generation. They parametrized the chemical feedback through a 
single quantity, Eg, which represents the kinetic energy input per unit gas 
mass of outflows from Pop III galaxies. This quantity is related to the num- 
ber of exploding Pop III stars and therefore encodes the dependence on the 
assumed IMF. For all values of the feedback parameter, Eg, Scannapieco 
et al. (2003) found that, while the peak of Pop III star formation occurs at 
z ^ 10, such stars continue to contribute appreciably to the star formation 
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rate density at much lower redshift, even though the mean IGM metal- 
Hcity has moved weU past the critical transition metallicity. This finding 
has important implications for the development of efficient strategies for 
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Figure 1.21. Top-heavy to normal IMF transition redshift, z/, as a function 
of SN-yT progenitor mass fraction and the mass density contributed by VMBHs, 
i^vMBH- Top: The computed critical density of VMBH remnants is compared 
to the observed values for SMBHs (upper dashed line). Bottom: The computed 
critical density of VMBH remnants is compared to the contribution to il from the 
X-ray-bright, off-center i^OSAT sources {lower dashed line) and to the abundance 
predicted assuming that the baryonic dark matter in galaxy halos is entirely con- 
tributed by VMBHs (upper dashed line). The observations on Qvmbh constrain 
the value of /^^. For a given /^^, the corresponding value for the transition 
redshift can be inferred by the Zf curve. This Figure is taken from Schneider et 
al. (2002). 



the detection of Pop III stars in primeval galaxies. At any given redshift, 
a fraction of the observed objects have a metal content that is low enough 
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to allow the preferential formation of Pop III stars. As metal-free stars are 
powerful Lya emitters (Tumlinson, Giroux & ShuU 2001; Schaerer 2003), 
it is natural to use this indicator as a first step in any search for primordial 
objects. Scannapieco et al. (2003) derived the probability that a given 
high-redshift Lya detection is due to a cluster of Pop III stars. In Figure 
ll.22l the isocontours in the Lya luminosity-redshift plane are shown, indi- 
cating this probability for various feedback efficiencies E^'^ . We see that 
Pop III objects populate a well-defined region, whose extent is governed by 
the feedback strength. Above the typical flux threshold, Lya emitters are 
potentially detectable at all redshifts beyond 5. Furthermore, the fraction 
of Pop III objects increases with redshift, independently of Eg'^ . For the 
fiducial case. Eg'' = 10"'^, the fraction is only a few percent at z = 4 but 
increases to approximately 15% by z = 6. So, the Lya emission from al- 
ready observed high-z sources can indeed be due to Pop III objects, if such 
stars were biased toward high masses. Hence collecting large data samples 
to increase the statistical leverage may be crucial for detecting the elusive 
first stars. 

1.6.3 Radiative Feedback and Reionization 

Before metals are produced, the primary molecule which acquires sufhcient 
abundance to affect the thermal state of the pristine cosmic gas is molecular 
hydrogen, H2. The main reaction able to produce molecular hydrogen are 
(Abel et al. 1997) 



H + c" ^ H" + hv, (1.54) 

H- -h H ^ H2 + e", (1.55) 

and 

H+ + H ^ Hj + hv, (1.56) 

H+ + H ^ H2 + H+, (1.57) 

In objects with baryonic masses > 3 x 10"* Mq, gravity dominates and 
results in the bottom-up hierarchy of structure formation characteristic of 
CDM cosmologies; at lower masses, gas pressure delays the collapse. The 
first objects to collapse are those at the mass scale that separates these two 
regimes. Such objects have a virial temperature of several hundred degrees 
and can fragment into stars only through molecular cooling (Tegmark et al. 
1997). In other words, there are two independent minimum mass thresholds 
for star formation: the Jeans mass (related to accretion) and the cooling 
mass. For the very first objects, the cooling threshold is somewhat higher 
and sets a lower limit on the halo mass of ~ 5 x 10"* Mq at z ~ 20. 
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Figure 1.22. Fraction of PopIII objects as a function of Lya luminosity and 
redshift. Isocontours of fractions > 10~ , 10~ '^,10"^ and 10 "^ are shown. 
Burst-mode star formation with a //^ = //^^ = 0.1 is assumed for all objects. In 
the PopIII case, a lower cutoff mass of 50Mq is assumed. Each panel is labeled 
by the assumed E^^ value. For reference, the dashed line gives the luminosity 
corresponding to an observed flux of 1.5 x lO"^'^ ergs cm~^ s~^, and the various 
points correspond to observed galaxies. The filled diamond is from Dey et al. 
(1998), the filled triangle is from Hu et al. (1999), the filled star is from Ellis et 
al. (2001), the open diamond is from Dawson et al. (2002), the open square is 
from Hu et al. (2002), the asterisks are from Lehnert & Bremer (2003), the open 
triangles are from Rhoads et al. (2003), the filled circles are from Fujita et al. 
(2003), and the filled squares are from Kodaira et al. (2003). The curves have 
been extended slightly past 2 = 4 for comparison with the Fujita et al. (2003) 
data-set. This Figure is taken from Scannapieco, Schneider & Ferrara (2003). 
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Figure 1.23. Equilibrium cooling curve at the center [solid lines) and virial 
radius [dashed line) of an isothermal halo at 2; = 9, as a function of virial tem- 
perature Tvir, and in the absence of a photoionizing background (i.e. prior to 
the reionization epoch). The halo has an assumed baryonic mass fraction of 
Qb = 0.019/i~^. The gas density dependence of the cooling function at high tem- 
peratures is due to Compton cooling off comic microwave background photons. 
The labeled curves extending to low temperatures show the contribution due to 
H2 for three assumed values of the metagalactic flux in the Lyman- Werner bands, 
47rJLw (in units of 10~^^ ergcm"^ s~^ Hz"^). This Figure is taken from Madau, 
Ferrara & Rees (2001). 
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Figure 1.24. The minimum mass needed to collapse. The function Mcizvir^ 
is plotted as a function on virialization redshift for standard CDM (SI = 1, 
Q.b = 0.06, h — 0.5). Only lumps whose parameters {zvir,M) lie above the 
shaded area can collapse and form luminous objects. The dashed straight lines 
corresponding to Tvir ~ 10* K and T^ir = 10"^ K are shown for comparison 
(dashed). The dark-shaded region is that in which no radiative cooling mecha- 
nism whatsoever could help collapse, since r„i,. would be lower than the CMB 
temperature. The solid line corresponds to 3 — cr peaks in standard CDM, normal- 



ized to (T8 = 0.7, so such objects with baryonic mass Sib x 2 x 10 Mq 
can form &t z = 30. This figure is taken from Tegmark et al. (1997). 
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As the first stars form, tlieir pliotons in tlie energy range 11.26-13.6 eV 
are able to penetrate the gas and photodissociate H2 molecules both in the 
IGM and in the nearest collapsing structures, if they can propagate that 
far from their source. Thus, the existence of an UV background below the 
Lyman limit due to Pop III objects, capable of dissociating the H2, could 
deeply influence subsequent small structure formation. 

Ciardi, Ferrara & Abel (2000) have shown that the UV flux from these 
objects results in a soft (Lyman - Werner band) UV background (SUVB), 
JlWj whose intensity (and hence radiative feedback efficiency) depends on 
redshift. At high redshift the radiative feedback can be induced also by the 
direct dissociating flux from a nearby object. In practice, two different sit- 
uations can occur: i) the collapsing object is outside the dissociated spheres 
produced by pre-existing objects: then its formation could be affected only 
by the SUVB {JLW,b), as by construction the direct flux {JLW,d) can only 
dissociate molecular hydrogen inside this region on time scales shorter than 
the Hubble time; ii) the collapsing object is located inside the dissociation 
sphere of a previously collapsed object: the actual dissociating flux in this 
case is essentially given by JLW.max = [JLW.b + JLW.d)- It is thus assumed 
that, given a forming Pop III, if the incident dissociating flux {JLW,b in 
the former case, JLW,max in the latter) is higher than the minimum flux 
required for negative feedback (Js), the collapse of the object is halted. 
This implies the existence of a population of "dark objects" which were 
not able to produce stars and, hence, light. 

To assess the minimum flux required at each redshift to drive the 
radiative feedback Ciardi et al. (2000) have performed non-equilibrium 
multifrequency radiative transfer calculations for a stellar spectrum (as- 
suming a metallicity Z = 10~^) impinging onto a homogeneous gas layer, 
and studied the evolution of the following nine species: H, H~, 11+ , He, 
Hc+, Hc++, H2, HJ and free electrons for a free fall time. We can then 
define the minimum total mass required for an object to self-shield from an 
external flux of intensity Js^ at the Lyman limit, as Mgh = (4/3) Tr{pii)R^j^, 
where (ph) is the mean dark matter density of the halo in which the gas 
collapses; Rsh is the shielding radius beyond which molecular hydrogen is 
not photodissociated and allowing the collapse to take place on a free-fall 
time scale. 

Values of Msh for different values of Jsfi have been obtained at various 
redshifts (see Figurc [T.25|l . Protogalaxies with masses above Mh for which 
cooling is predominantly contributed by Lya line are not affected by the 
radiative feedback studied here. The collapse of very small objects with 
mass < Merit is on the other hand not possible since the cooling time is 
longer than the Hubble time. Thus radiative feedback is important in the 
mass range 10^ — IO^Mq, depending on redshift. In order for the negative 
feedback to be effective, fluxes of the order of 10~^^ — 10~^^erg s~^ cm~^ 
Hz^^ sr~^ are required. These fluxes are typically produced by a Pop III 
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with baryonic mass lO^Af^ sMq at distances closer than ~ 21 — 7x Af^ g kpc 
for the two above flux values, respectively, while the SUVB can reach an 
intensity in the above range only after z ra 15. This suggests that at 
high z negative feedback is driven primarily by the direct irradiation from 
neighbor objects in regions of intense clustering, while only for z < 15 the 
SUVB becomes dominant. 
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Figure 1.25. Minimum total mass for self-shielding from an external incident 
flux with intensity Js,o at the Lyman limit. The curves are for different redshift: 
from the top to the bottom z =8.3, 11.4, 15.4, 18.0, 22.1, 29.6. Circles show the 
value of Mh (open) and Merit (filled). Radiative feedback works only in the solid 
portions of the curves. This Figure is taken from Ciardi et al. (2000) 



Figure [T.26l illustrates all possible evolutionary tracks and final fates of 
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primordial objects, together with the mass scales deterrained by the various 
physical processes and feedbacks. There are four critical mass scales in the 
problem: (i) Merit, the minimum mass for an object to be able to cool in a 
Hubble time; (ii) Mh, the critical mass for which hydrogen Lya line cooling 
is dominant; (iii) Msh, the characteristic mass above which the object 
is self-shielded, and (iv) Mty the characteristic mass for stellar feedback, 
below which blowaway can not be avoided. Starting from a virialized dark 
matter halo, condition (i) produces the first branching, and objects failing 
to satisfy it will not collapse and form only a negligible amount of stars. In 
the following, we will refer to these objects as dark objects. Protogalaxies 
with masses in the range Merit < M < Mh are then subject to the effect of 
radiative feedback, which could either impede the collapse of those of them 
with mass M < Mgh, thus contributing to the class of dark objects, or allow 
the collapse of the remaining ones (M > Msh) to join those with M > Mh 
in the class of luminous objects. This is the class of objects that convert 
a considerable fraction of their baryons in stars. Stellar feedback causes 
the final bifurcation by inducing a blowaway of the baryons contained in 
luminous objects with mass M < M^y-, this separates the class in two 
subclasses, namely "normal" galaxies (although of masses comparable to 
present day dwarfs) that we dub gaseous galaxies and tiny stellar aggregates 
with negligible traces (if any) of gas to which consequently we will refer to 
as naked .stellar clusters. 

The relative numbers of these objects are shown in Figure IT. 271 The 
straight lines represent, from the top to the bottom, the number of dark 
matter halos, dark objects, naked stellar clusters and gaseous galaxies, 
respectively. The dotted curve represents the number of luminous objects 
with large enough mass (M > Mh) to make the H line cooling efhcient and 
become insensitive to the negative feedback. We remind that the naked 
stellar clusters are the luminous objects with M < M^y, while the gaseous 
galaxies are the ones with M > Mby] thus, the number of luminous objects 
present at a certain redshift is given by the sum of naked stellar clusters 
and gaseous galaxies. We first notice that the majority of the luminous 
objects that are able to form at high redshift will experience blowaway, 
becoming naked stellar clusters, while only a minor fraction, and only at 
z <;15, when larger objects start to form, will survive and become gaseous 
galaxies. An always increasing number of luminous objects is forming with 
decreasing redshift, until z ~ 15, where a flattening is seen. This is due 
to the fact that the dark matter halo mass function is still dominated 
by small mass objects, but a large fraction of them cannot form due to 
the following combined effects: i) toward lower redshift the critical mass 
for the collapse (Merit) increases and fewer objects satisfy the condition 
M > Merit', ii) the radiative feedback due to either the direct dissociating 
flux or the SUVB increases at low redshift as the SUVB intensity reaches 
values significant for the negative feedback effect. When the number of 
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Figure 1.26. Possible evolutionary tracks of objects as determined by the pro- 
cesses and feedbacks included in the model. See text. This figure is taken from 
Ciardi et al. (2000) 
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luminous objects becomes dominated by objects with M > Mjj, by z ^^ 10 
the population of luminous objects grows again, basically because their 
formation is now unaffected by negative radiative feedback. A steadily 
increasing number of objects is prevented from forming stars and remains 
dark; this population is about ^ 99% of the total population of dark matter 
halos at z ~ 8. This is also due to the combined effect of points i) and ii) 
mentioned above. This population of halos which have failed to produce 
stars could be identified with the low mass tail distribution of the dark 
galaxies that reveal their presence through gravitational lensing of quasars. 
It has been argued that this population of dark galaxies outnumbers normal 
galaxies by a substantial amount, and Figure H . 271 supports this view. 

1.6.3.1 Escape of Ionizing Photons 

One crucial point to understand cosmic reionization is to determine the 
fraction of ionizing photons that can escape from their parent galaxy. By 
solving the time-dependent radiation transfer problem of stellar radiation 
through evolving superbubbles within a smoothly varying H I distribution, 
Dove et al. (2000) have estimated the fraction of ionizing photons emit- 
ted by OB associations that escapes the H I disk of our Galaxy into the 
halo and intergalactic medium (IGM). They considered both coeval star- 
formation and a Gaussian star-formation history with a time spread at — 2 
Myr; the calculations are performed both a uniform H I distribution and a 
two-phase (cloud/intercloud) model, with a negligible filling factor of hot 
gas. The shells of the expanding superbubbles quickly trap or attenuate 
the ionizing flux, so that most of the escaping radiation escapes shortly 
after the formation of the superbubble. Superbubbles of large associations 
can blow out of the H I disk and form dynamic chimneys, which allow 
the ionizing radiation to escape the H I disk directly. However, blowout 
occurs when the ionizing photon luminosity has dropped well below the 
association's maximum luminosity. For the coeval star-formation history, 
the total fraction of Lyman Continuum photons that escape both sides of 
the disk in the solar vicinity is {fesc) — 0.15 ± 0.05 (Figure ri.28|) . For the 
Gaussian star formation history, {fesc) — 0.06 ± 0.03. 

Bianchi, Cristiani & Kim (2001) derived the Hl-ionizing background, 
resulting from the integrated contribution of QSOs and galaxies, taking 
into account the opacity of the intervening IGM. They have modeled the 
IGM with pure-absorbing clouds, with a distribution in column density of 
neutral hydrogen, TVhi, and redshift, z, derived from recent observations of 
the Lya forest (Kim et al. 2002) and from Lyman Limit systems. The QSOs 
emissivity has been derived from the recent fits of Boyle et al. (2000), while 
the stellar population synthesis model of Bruzual & Chariot (2003) and a 
star-formation history from UV observations for the galaxy emissivity. Due 
to the present uncertainties in models and observations, three values for the 
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Figure 1.27. Number evolution of different objects in the simulation box. This 
Figure is taken from Ciardi et al. (2000) 



fraction of ionizing photons that can escape a galaxy interstellar medium, 
fesc= 0.05, 0.1 and 0.4, as suggested by local and high-z UV observations 
of galaxies, respectively, are used. 

In Figure [T. 291 is shown the modeled UV background, J{i',z), at the 
Lyman limit as a function of redshift (solid lines) for the flat universe with 
r2m=l. The total background is shown as the sum of the QSO contribution 
(the same in each model; dotted line) and the galaxy contribution (scaled 
with fesc] dashed lines). 

At high redshift, the value of the UV background is constrained by 
the analysis of the proximity effect, i.e. the decrease in the number of 
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Figure 1.28. Photon luminosity emitted by a single OB association (solid line) 
and the luminosity of photons escaping each side of the H I disk (dashed line) for 
coeval star-formation. Also shown is the fraction of photons emitted that escape, 
VesciNT,t) = Sesc{t)/S{t). Here, Ztr = 2(7^, the cavity height is Zcy = Agu, 
and 13=1. The I-front is trapped within the shell during the time interval 
where rjesc = 0. After blowout, rjesc suddenly increases even though the photon 
luminosity of the association is relatively low at these times. Case with A''* = 200. 
This Figure is taken from Dove et al. (2000) 
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intervening absorption lines that is observed in a QSO spectrum when 
approaching the QSO's redshift (Bajthk et aL 1988). Using high res- 
olution spectra, Giallongo et al. (1996) derived J(912A) = 5.0l^-^ x 
10~^^ erg cm^^ s^^ Hz^^ sr~^ for 1.7 < z < 4.1 (see also Giallongo et 
al. 1999). Larger values are obtained (same units) by Cooke et al. (1997), 
J(912A) = l.OtoJ ^ 10"^^ for 2.0 < z < 4.5 and by a recent re-analysis of 
moderate resolution spectra by Scott et al. (2000) who found J(912A) = 
7.0^4^4 -10^^^ in the same redshift range. These measurements are shown in 
Figure 11.291 with a shaded area: the spread of the measurements obtained 
with different methods and data gives an idea of the uncertainties associ- 
ated with the study of the proximity effect. Several models can produce a 
value of J(912A) compatible with one of the measurements from the prox- 
imity effect at z '^ 3 shown in Figure 11.291 from a simple QSO-dominated 
background to models with fesc ^ 0.2, leading to a limit of fesc ^ 20%. 

1.6.3.2 Cosmological Ionization Fronts 

The radiative transfer equation in cosmology describes the time evolution 
of the specific intensity I^, = I{t, x, uj, v) of a diffuse radiation field (Ciardi, 
Ferrara& Abel 2000): 



Idh n- V/^ H{t) 

C dt Oem C 



V— 31^ 



V.-Xh, (1.58) 



where c is the speed of light, x is the continuum absorption coefficient per 
unit length along the line of sight, and 77,^ is the proper space averaged 
volume emissivity. 

If massive stars form in Pop III objects, their photons with hiy > 
13.6 eV create a cosmological HII region in the surrounding IGM. Its radius, 
Ri, can be estimated by solving the following standard equation for the 
evolution of the ionization front: 



^^'^ HR,^ ' 



dt ' i-KUHRf 



S^iO) - ^nR^nj^a^'^ 



(1.59) 



note that ionization equilibrium is implicitly assumed. H is the Hubble 
constant, 5*^(0) is the ionizing photon rate, nn is the IGM hydrogen number 
density and a^^-* is the hydrogen recombination rate to levels > 2. For 
aR ^ c/H the cosmological expansion term HRi can be safely neglected. 
In its full form equation (jl.59|l must be solved numerically. However, if 
steady-state is assumed {dRi/dt ~ 0), then Ri is approximately equal to 
the Stromgren radius R^z) = [{MN^/dt)/{A'Ka2n^H\z))Y^^: 

R,c^R,^ 0.05(r2f,/i2)-2/3(i + z)-^Sl(\pc, (1.60) 
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Figure 1.29. UV background at A = 912A for the models with /esc=0.05, 0.1 
and 0.4 (solid lines), in a flat flm ~ 1 universe. Also shown are the separate 
contribution of the QSOs (dotted line) and of the galaxies (dashed lines, each 
corresponding to a value of fesc)- The shaded area refer to the Lyman limit UV 
background estimated from the proximity effect (Giallongo ot al. 1996; Cooke 
et al. 1997; Scott et al. 2000). The arrow shows an upper limit for the local 
ionizing background (Vogel et al. 1995). The datapoint at z = 3 is derived from 
a composite spectrum of Lyman-break galaxies (Steidel et al. 2001). The models 
and the datapoint of Steidel et al. (2001) have been multiplied by a z-dependent 
factor, to take into account the cloud emission. This Figure is taken from Bianchi 
et al. (2001) 



where S'47 = S'i(0)/10^^s^^. The approximate expression for Si(0) is ~ 
1.2 X 10*^ fbMe s~^ when a 20% escape fraction for ionizing photons is 
assumed and fb is the fraction of baryon that is able to cool and become 
available to form stars. In general Rs represents an upper limit for Ri, since 
the ionization front completely fills the time- varying Stromgren radius only 
at very high redshift, z ~ 100. 
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1.6.3.3 Reionization after WMAP 

The WMAP satellite has detected an excess in the CMB TE cross-power 
spectrum on large angular scales {I < 7) indicating an optical depth to 
the CMB last scattering surface of t^ = 0.17 ± 0.04 (Kogut et al. 2003). 
The reionization of the universe must therefore have begun at relatively 
high redshift. Ciardi, Ferrara & White (2003) have studied the reioniza- 
tion process using supercomputer simulations of a large and representative 
region of a universe which has cosmological parameters consistent with the 
WMAP results. The simulation follow both the radiative transfer of ion- 
izing photons and the formation and evolution of the galaxy population 
which produces them. 

They have shown that the WMAP measured optical depth to electron 
scattering is easily reproduced by a model in which reionization is caused 
by the first stars in galaxies with total masses of a few x 10^ Mq. Moreover, 
the first stars are 'normal objects', i.e. their mass is in the range 1-50 Mq, 
but metal-free. Among the different model explored, the 'best' WMAP 
value for Te is matched assuming a moderately top-heavy IMF (Larson 
IMF with Mc = 5 Mq) and an escape fraction of 20%. A Salpeter IMF 
with the same escape fraction gives Te = 0.132, which is still within the 
WMAP 68% confidence range. Decreasing fesc to 5% gives Te — 0.104, 
which disagrees with WMAP only at the 1.0 to 1.5a level. In Figure IT. 301 
is shown the neutral hydrogen number density at three redshift {z = 17.6, 
15.5, and 13.7) and in Figure IT. 311 is shown the evolution of the electron 
optical depth for the different IMFs considered. 

In the best-fit model reionization is essentially complete by Zr ~ 13. 
This is difficult to reconcile with observation of the Gunn-Peterson effect in 
z > 6 quasars (Becker et al. 2001; Fan et al. 2002) which imply a volume- 
averaged neutral fraction > 10"'^ and a mass-averaged neutral fraction 
~ 1% at z = 6. Then a fashinating (although speculative) possibility is 
that the universe was reionized twice (Gen 2003; Whyte & Loeb 2002) with 
a relatively short redshift interval in which the IGM became neutral again. 
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Figure 1.30. Slice through the simulation boxes. The six panel show the neutral 
hydrogen number density for the case of a Larson IMF with Mc — 5 Mq (upper 
panels) and for a Salpeter IMF (lower panels) at redshifts, from left to right, 
z = 17.6, 15.5 and 13.7. The escape fraction is 20%. The box for the radiative 
transfer simulation has a comoving length of L = 20h~^ Mpc. This Figure is 
taken from Ciardi et al. (2003). 
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Figure 1.31. Redshift evolution of the electron optical depth, Te, for an Salpeter 
IMF with /esc = 5% (long-dashed line) and fesc ~ 20% (short-dashed line), and 
for a Larson IMF with Mc = 5 M© and /e^c = 20%. The dotted line refers 
to a sudden reionization aX z = 16. The shaded region indicates the optical 
depth Te = 0.16 ± 0.04 (68% CL) implied by the Kogut et al. (2003) 'model 
independent' analysis. In the inset the redshift evolution of the volume-averaged 
ionization fraction, x.^, is shown for the three runs. This Figure is taken from 
Ciardiet al. (2003). 
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